One- and two-atom states in a rotating ring lattice 
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We study the states of one and two atoms in a rotating ring lattice in a Hubbard type tight- 
binding model. The model is developed carefully from basic principles in order to properly identify 
the physical observables. The one-particle ground state may be degenerate and represent a finite 
flow velocity depending on the parity of the number of lattice sites, the sign of the tunneling matrix 
element, and the rotation speed of the lattice. Variation of the rotation speed may be used to control 
one-atom states, and leads to peculiar behaviors such as wildly different phase and group velocities 
for an atom. Adiabatic variation of the rotation speed of the lattice may also be used to control the 
state of a two-atom lattice dimer. For instance, at a suitably chosen rotation speed both atoms are 
confined to the same lattice site. 

PACS numbers: 03.75.Lm, 37.10.Jk, 05.30.Jp, 05.50.+q 
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I. INTRODUCTION 

Optical lattices have enhanced both AMO and con- 
densed matter physics, and will undoubtedly continue to 
do so for a long time. Theoretical analyses of optical lat- 
tices routinely resort to periodic boundary conditions, as 
if the lattice folded back onto itself into a ring. Usually 
this is just a matter of convenience, but the topology of 
a torus may have a profound effect on the physics. As 
an example, the phase winding of a supcrfluid around 
a ring cannot change discontinuously, which ultimately 
stabilizes a persistent current. Ring traps for atoms have 
in fact been demonstrated [1-4] and employed to study 
supcrfluid flow [1, 3]. An intriguing technique whereby 
a rapidly moving optical trap, when averaged over time, 
may "paint" not only a ring trap but also a virtually arbi- 
trary time dependent structure, in particular, a rotating 
ring lattice, has also been demonstrated [5]. A ring lat- 
tice with a precisely set and possibly small number of 
sites rotating at a controllable speed is within the reach 
of current experimental techniques. 

The other background clement here is the lattice dimcr 
made of two atoms in an optical lattice. Both the site- 
to-site tunneling and the atom-atom interactions could 
be controlled, by adjusting lattice parameters and with 
the aid of a Feshbach resonance. A lattice dimer may 
thus make a tailored custom molecule. The experimen- 
tal demonstration of a dimer bound by repulsive atom- 
atom interactions [6] has in part motivated theoretical 
work from several groups [6-16]. As usual, we have ap- 
plied periodic boundary conditions in our analysis [14- 
16], which is convenient and permissible when the lattice 
is long and the boundary conditions cannot matter in 
practice. However, the boundary conditions are impor- 
tant if the lattice truly is a finite-size ring, and some of 
our technical assumptions such as the even number of lat- 
tice sites [14] need to be re-examined. Besides, rotation 
of the lattice could provide another handle for controlling 
the molecules. 



In the present paper we first, Sec. II, study systemat- 
ically the effects of the rotation of the lattice on a Hub- 
bard type mode by expanding on our earlier coordinate- 
transformation arguments [17]. This groundwork allows 
us to identify the physical observable for one and two 
atoms in a rotating ring lattice, Sees. Ill and IV. Ther- 
mal preparation and adiabatic variation of the rotation 
speed prove to be effective methods to control the state 
of both one- and two-atom systems, and lead to quite a 
few perhaps surprising results. For instance, by slowly 
spinning up the lattice one cannot change the speed with 
which an atom emerges after it is released from the lat- 
tice, but a localized atomic wave packet will almost track 
the varying rotation speed; and at certain rotation speeds 
a bound dimer of two atoms is confined to a single lattice 
site. We conclude in Sec. V with a few brief comments. 



II. HUBBARD TYPE MODELS 

Several methods have been used in the past to address 
the effects of rotation on optical lattices [17-19], and by 
now the observation that the rotation leads to phase fac- 
tors in site-to-site tunneling matrix elements has the force 
of folklore. We find the same basic result, but emphasize 
that one needs to keep track of the physical observables 
of the system carefully. 

One may resort to the formal similarity between ro- 
tation and magnetic field, say, in that the Coriolis force 
derives from a vector potential, and argue basically in 
terms of minimal-coupling substitution [18]; develop for 
the lattice an approximation to the usual — o;-L term that 
emerges in a transformation to the rotating frame [19]; 
or analyze coordinate transformations directly [17, 20]. 
All of these methods have their problems, however. For 
instance, in quantum mechanics coordinate transforma- 
tions are mathematically unitary transformations and 
change the appearance of quantum mechanics. Similarly, 
minimal substitution is carried out in a given fixed gauge, 



but in quantum mechanics a transformation to a differ- 
ent gauge is also a unitary transformation. The following 
question exemplifies the core of the problem: In which 
unitarily transformed version of quantum mechanics is 
e jpr the wave function representing a particle with the 
velocity p/m? The debates about the p-A and d-E forms 
of the dipole interaction that have erupted periodically 
in the past are a manifestation of this type of ambigu- 
ity [21]. We resolve such issues by carefully expanding on 
our earlier coordinate-transformation approach [17]. In 
the end, we will keep track of three reference frames. 

We start from the assumption that the motion of the 
atoms is confined to a torus with the circumference L. In 
the limit of an asymptotically strong transverse confine- 
ment, the coordinate along the direction of the torus, x, 
remains the only relevant degree of freedom. In principle 
one might consider angular momenta, but for sufficiently 
tight transverse confinement we may model the motion 
in the x direction simply as one-dimensional translation. 
The physics takes place over the interval [— L/2, L/2), 
but here we imagine that the coordinate x ranges over 
the entire real axis and impose periodic boundary condi- 
tions over the distance L to account for the topology of 
the torus. The inner product of one-atom wave functions 
is defined as 



(iM) 



L/2 



dxip*(x)(f>(x) . 



(1) 



-L/2 



from which the Hilbert space structure follows. 

Suppose next that we have a stationary potential V(x) 
along the direction of the torus. We take the potential 
V(x) to be periodic over the distance L and also over 
a shorter distance a = L/N, so that we have an TV-site 
lattice with the lattice spacing a and periodic boundary 
conditions from end to end of the lattice. We denote the 
sites of the lattice, e.g., minima of the potential V(x), by 
x n , and for future reference also define the wave vectors 
fc„. We set 



A',i 



2-KTl 



(2) 



The choice here (by no means unique) is that, unless 
otherwise specified, for an even number of the sites in 
the lattice the index n ranges from —N/2 to N/2 — 1 in 
unit steps, while for an odd number of lattice sites N 
the range is from -N/2 - 1/2 to N/2 - 1/2. The wave 
vectors thus run over the first Brillouin zone of solid-state 
physics. 

Assume now that the potential V(x) is made to rotate 
along the torus so that we have a rotating ring lattice, 
with the velocity along the ring specified by v. The way 
that we proceeded in Rcf. [17] is probably uncontrover- 
sial, but a closer examination reveals a number of sub- 
tleties. 

To begin with, the notion of an explicitly time depen- 
dent potential energy is not part of the standard edifice 
of classical mechanics, and so not part of quantum me- 
chanics either. We resolve this immediate issue with the 



Galilean invariance principle that in a frame moving with 
the potential energy at the uniform speed v, the classi- 
cal one-particle physics is the same as in the stationary 
frame when the potential does not move. In terms of the 
usual classical position and conjugate momentum vari- 
ables, the Hamiltonian is 



H v = ^- + V(x) . 
2vn 



(3) 



We may quantize the Hamiltonian in the usual way by 
replacing the position and momentum variables with the 
quantum operators x and p, to have 



H„ 



P 

2m 



V(x). 



(4) 



Interestingly, the canonical commutator does not 
uniquely determine the quantized momentum operator 
p. For instance, if [x,p] = ih, then [x,p+ f(x,t)] = iH 
also holds true for an arbitrary function /(x,i). For our 
purposes it suffices to set f{x,t) = —A, a so far undeter- 
mined constant. In position representation, in the frame 
moving with the lattice, we therefore write the quantum 
Hamiltonian as 



2m \i dx 



+ V[x) + K 



(5) 



Here we have added another constant K that has no ef- 
fect on the dynamics, but is at our disposal for later 
convenience. 

Let us now transform from the moving coordinate sys- 
tem to the stationary frame. The coordinate transforma- 
tion, a Galilean transformation, is defined by 



t. 



£■ 



(6) 



where £ and r are position and time in the stationary lab 
frame. It follows from the time dependent Schrodinger 
equation in the rotating frame for the wave function 
ijj v (x,t) that the corresponding stationary-frame wave 
function 



V>(£,t) = i() v (£ + vt,t) 



(7) 



also satisfies the time dependent Schrodinger equation 
with the Hamiltonian 



H = ± 
2m 



h d (X 



+ V(£ - vt) 



+v\ — i mi) 2 



K . 



(8) 



The crux of the present argument is as follows: In the 
stationary frame, and in the absence of the moving poten- 
tial (V — > 0), we want to regain the ordinary Schrodinger 
equation for a free particle, with the usual free-particle 
Hamiltonian 



H 



W_cP_ 

2~^d^ 



(9) 



and periodic boundary condition for the wave function 
?/>(£, t) in the variable £. This happens only if we choose 



X = mv, K = -^mv z . (10) 

This puts the rotating- frame Hamiltonian into the form 
1 (h d 



2m \ i dx 



V(x) 



\ mv 



(11) 



The rotating- frame wave function ij; v (x,t) has periodic 
boundary conditions in the variable x. 

The Hamiltonian (11) is still nonstandard form. To 
remedy this situation, we introduce a transformation of 
the wave function ip v — > \IJ, a momentum translation, by 



i> v (x,t) 



*0M) 



(12) 



Upon this unitary transformation the Hamiltonian be- 
comes 

h 2 a 2 



Hi 



2m dx 2 



+ V{x)-\mv 2 , (13) 



and the wave function ^> acquires twisted boundary con- 
ditions, 



$(x + L,t) = e^ %m Tr- qj( x , t) 



(14) 



Our development has introduced three frames: The 
stationary lab frame (LF), the rotating frame (RF), and 
the momentum translated rotating frame (MTRF). In 
quantum mechanics the frames are related by unitary 
transformations. By choice, the bulk of our analysis is in 
the MTRF, but all three frames figure in the argument 
and have to be distinguished carefully. 

We illustrate the distinction between these frames by 
solving for the stationary states in the absence of any 
external potential, V = 0, in the MTRF. The process 
is standard. In particular, the operator 5 -^ is hermi- 
tian with respect to the inner product (l) even with 
the MTRF twisted boundary conditions. The normal- 
ized stationary states, solutions to the time dependent 
Schrodinger equation for which all quantum mechanical 
expectation values are independent of time, read 



* n (x,t) = -^=e< K » 



-n„t) 



(15) 



with 



K — h — h • 

k v = mv/h; k r , 



27m 



n = 0,±l,.. 



(16) 



To go from the MTRF to the LF we apply both the 
appropriate momentum translation and the inverse of 
the Galilean transformation (6), and find the stationary 
states 



ifci") 



1 



VI 



e i(fc„e-a; n T) 



hkl 
2m 



(17) 



These are the usual stationary states of the LF Hamil- 
tonian (9), and satisfy the proper LF periodic boundary 
conditions. 

Simple as this exercise is, it already underscores the 
problem that one encounters in connection with statis- 
tical mechanics. Namely, the energies of the states in 
the MTRF and in the LF arc different, Ml n and Huj n , 
and the difference is not just a constant. What would be 
the thermal equilibrium, and even the ground state that 
gets prepared at zero temperature, is therefore ambigu- 
ous. The question becomes even more acute if there is a 
potential present, because the Hamiltonian then depends 
explicitly on time in the stationary frame and standard 
statistical mechanics is inapplicable. We will not attempt 
to clean up these issues, but simply hypothesize a reso- 
lution as we go along. 

We next proceed to Hubbard type models for a lattice. 
Here and in the rest of the paper we restrict ourselves to 
the basic one-band (tight-binding) model, and ignore al- 
together the possible coupling of the energy bands that 
may occur with increasing strength of atom-atom inter- 
actions [22, 23]. 

We begin in the MTRF, Hamiltonian (13), with what 
is essentially a rederivation of Bloch's theorem. Thus, let 
us introduce the lattice translation operator U and its 
inverse U~ l via their actions on the wave function by 

{U^){x)=^(x + a), (U- 1 y)(x) = y(x-a). (18) 

This operator preserves the inner products, (t/\P, U$) = 
(W, $), and is invertiblc, so it is unitary. Moreover, for a 
periodic potential V(x) it commutes with the Hamilto- 
nian H' v of Eq. (13). Therefore H' v and U can be diago- 
nalized simultaneously. Let \P be a simultaneous eigen- 
function and A the corresponding eigenvalue of U . By 
virtue of the twisted boundary conditions we have 

(U N ^)(x) = \ N V(x) = V(x + L) = e- tk " L y(x) . (19) 

The possible eigenvalues of U are therefore of the form 



A, 



iK n a 



(20) 



see Eqs. (16) for the notations used in the present argu- 
ment. By writing an eigenvalue and eigenvector corre- 
sponding to K n in the form 



*„(x) = e lK " x u n (x) 



(21) 



we immediately sec that u n (x) has to be periodic over the 
distance a. This is Bloch's theorem for twisted boundary 
conditions. 

Obviously, by virtue of the unitary transforma- 
tion (12), we would write Bloch's theorem in the RF, 
for the Hamiltonian (11), in the form 

jp n (x) = e^ x u n (x) , (22) 

with the same functions u n {x) as in Eq. (21). 



When the lattice potential is weak, the rotation may 
rearrange the band structure. However, as we are in- 
terested in the tight-binding limit, we assume that the 
lowest energy band remains essentially undeformed by 
rotation, and consider only the N states therein. 

Let us next define the analogs of Wannier functions, 
given the twisted boundary conditions: 

v n' 

= -L e~ iKx Y^ elK ' (x ~ Xn) u n > (x) 

J_ e -ik v x J- e^' (*-*»W (x - x n ) 

(23) 
(24) 



N 



AkvX w{x - X n ) 



-ik v x r 



W(x-x n ), 



where 



(x) = -^=^2e lk " x u n (x) 



(25) 



is evidently the Wannier function in the RF for the lattice 
position n = 0, and 



W(x) = e- iKx w(x) 



(26) 



likewise in the MTRF. Incidentally, an attempt to con- 
struct Wannier functions numerically quickly reveals that 
they are not unique. Namely, the Bloch functions (21) 
can be multiplied by arbitrary phase factors to give 
equally good Bloch functions, but possibly completely 
different Wannier functions. Our usual choice is to put 
a minimum of the lattice potential at x = and pick 
the Bloch states so that \I / n (0) is real and positive. Em- 
pirically, this choice seems to produce Wannier functions 
W n (x) , each of which is narrowly centered at x n . What- 
ever the choice, according to Eq. (24) the Wannier func- 
tions W n {x) are translated copies of a single function 
W(x). The functions W n (x) are orthonormal if the orig- 
inal Bloch states (21) are, and in fact make an orthonor- 
mal basis for the states in the lowest energy band. 

The MTRF Wannier functions W n (x) and W(x) = 
Wo(x) themselves obey twisted boundary conditions by 
construction, and so do all of their linear combina- 
tions. However, Eq. (24) also suggests another mecha- 
nism whereby twisted boundary conditions may be sat- 
isfied. Namely, a linear combination of functions of this 
form with an infinite number of coefficients c„ satisfy- 
ing periodic boundary conditions, c„ = c„+at, if it ex- 
ists in the first place, satisfies twisted boundary condi- 
tions no matter what the function W(x) is. Conversely, 
if some W(x) were strictly localized around x = 0, equal 
to zero by |a;| > L/2, and did not repeat at the inter- 
vals of L, then the only way to make a wave packet sat- 
isfying twisted boundary conditions from the functions 



W n {x) 



-ik v x r 



W(x — x n ) would be to insist on peri- 



odic boundary conditions for the expansion coefficients 



In short, the following modeling suggests itself: 

(i) We adopt a function W(x) that is localized around 
x = 0, without periodic recurrences, and in such a 



way that W n (x) 



-ik v x r 



W{x — x n ) for different 



n may be reasonably regarded as an orthonormal 
set of functions. 

(ii) We consider only linear combinations of the func- 
tions W n (x) of the form 



*(a:)= Yl c nW n (x) 



E 



c n c 



-ik v x n 



W(x 



(27) 



where the expansion coefficients satisfy periodic 
boundary conditions, c n = c n +N- 

(iii) The wave function W n (x) = e~ lk '° Xn W [x — x n ) is 
the one-particle state representing an atom that re- 
sides at the site n, and b n is the corresponding bo- 
son operator. 

The function W(x) specifying the Wannier functions 
restricted to the interval [—L/2, L/2) provides an exam- 
ple of a function of this kind. Specifically, given the 
Wannier function around x — 0, Wq(x), an d the unit 
step function 9{x), we could define W{x) = Wq(x)6(x + 
L/2)9(L/2 — x). The ensuing construction of the lattice 
states is then precisely the same as if we simply consid- 
ered the original Wannier function restricted to a proper 
finite set of indices such as n € [— N/2, N/2 — 1). Like- 
wise, one could use an approximation to the ground state 
of a particle in the potential well at x = to model the 
function W{x). Such constructions are most useful if the 
dependence on rotation is primarily contained in the pref- 
actors e ~ tk,,Xn of the basis wave functions W n (x). After 
we have adopted these one-particle states and boson op- 
erators, the lattice per se is periodic, with the sites n and 
n + N regarded as the same. 

Given the functions W n {x), we finally set up the corre- 
sponding Hubbard model with nearest-neighbor tunnel- 
ing and atom-atom interactions in the usual manner, 

f = \ E H e# Jbi^bn+e-^rbiA^+ubib^A. 

n 

(28) 
Logic would dictate that we use the notation H' v for this 
MTRF operator, but we write H anyway. We have 

2 r L/2 

J = — -/ dxW*{x-a){H' v W){x) (29) 

" J -L/2 

2 r L/2 

-e-* / dxw*(x-a){H v w)(x). (30) 



h 



L/2 



Since the periodic boundary conditions of the MTRF 
are awkward to deal with, using the transformation (12) 



we have also expressed the tunneling matrix element in 
terms of what would approximate the Wannier functions 
for the RF Hamiltonian (11). The phase factor <j> equals 
the phase twist per lattice site, or 






$ 



mvL 



(31) 



where $ is the end-to-end phase twist over the lattice. 
The atom-atom interaction part depends on the s-wave 
scattering length ao, but to obtain it accurately we need 
to know also about the structure of the wave function in 
the directions transverse to x. Given a full 3D form of 
the Wannier functions, we could write 



V 



Airhari 



d 3 x\W(x)\ 4 . 



(32) 



A periodic lattice is implied, so that the sites n = 
and n = N are the same. We have finally dropped the 
constant — ^mv 2 in the Hamiltonian as it has no effect 
on either dynamics or thermodynamics. 

We also need the Hamiltonian in lattice momentum 
representation. Specifically, let us define the boson oper- 
ators in lattice momentum space B q in such a way that 

& » = 7f E eiqnB « * B « = 7^ E e_i9 " 6 « > ( 33 ) 



then the Hamiltonian reads 



H = -jJ2cos(q-<f>)BlB q 



U 
"2N 



E S qi+q2,q3+q4 B qi B q 2 B<l3 B q4 ■ ( 34 ) 

91,92,93,94 

With definitions in terms of the lattice spacing a such as 
q = ak n , the discrete numbering of the lattice momenta 
is hidden and the lattice momenta q are made dimen- 
sionlcss. Moreover, by virtue of the periodicity in lattice 
momentum space, sums and comparisons of lattice mo- 
menta are modulo 2tt unless there is an explicit reason to 
proceed differently. If the original operators b n annihilate 
Wannier states of the lattice, the operators B q annihilate 
Bloch states. 

Several items remain to be cleaned up. First, some 
fine tuning on the rotation phase is in order. Suppose 
that we have at the origin the harmonic oscillator poten- 
tial V(x) = imw 2 i 2 , and use the ground state following 
from the RF Hamiltonian H v , Eq. (11) as the "Wannier 
function" w(x) in Eq. (30). The tunneling matrix ele- 
ment is then 



J 



we 



(35) 



independent of the rotation phase <f>. On the other hand, 
when we use a sinusoidal potential V(x) as appropri- 
ate for an ideal optical lattice in a numerical solution of 
the RF Hamiltonian (11), we find that both the magni- 
tude and the phase of the tunneling amplitude J depend 



somewhat on cj>. In particular, there is a phase drag of 
sorts incorporated into J. However, in the tight-binding 
limit these are small corrections. Given that our model 
is heavy-handed in many other details anyway, we will 
assume that the phase <f> accounts for all effects of the 
rotation of the lattice. 

Second, one could augment each of the functions 
W n (x) with an arbitrary phase factor e~ tiPn , 



W n (x) 



>W n (x) 



(36) 



Indeed, this operation does not change the orthonormal- 
ity of the functions W n (x), or the function space they 
span. The end result would be additional phase factors 
in the site-to-site tunneling amplitudes, such as 



e*+JbL,h 



n+l u n 



ju n+l u n ■ 



(37) 



The lattice remains periodic, wraps around at n = N, 
only if the added phase factors wrap around, so that we 
have e' llPn+N = e lVn . This means that the phases added to 
the tunneling matrix elements, <p n +i ~ <Pm must add up 
to an integer multiple of 2tt. Since the absolute phases 
assigned to the members of an orthonormal basis can- 
not have any effect on the physics, the total transfor- 
mation (36) and (37) to alter the phase winding cannot 
change the physics of the lattice in any way. 

Because the effect of the transformation of Eqs. (36) 
and (37) on the Hamiltonian (28) with the choice ip n — 
nA<fr/N is exactly the same as the change of the rotation 
phase by A</>, one might be tempted to conversely sur- 
mise that the physics repeats periodically when the rota- 
tion speed is increased; is unchanged whenever $ changes 
by 2ir, or the rotation phase <f> changes by 2ir/N. This 
notion, if valid, would severely limit the possibilities to 
control the lattice by varying the rotation speed. Fortu- 
nately, it is incorrect. 

The transformation (36) comes with both phase 
changes in the basis functions and an apparent change in 
the rotation phase (37) that cancel inasmuch as physical 
observables are concerned. On the other hand, a change 
in the rotation speed and the accompanying changes in 
the Hamiltonian do not per se alter the basis functions. 
The spectrum of any observable, such as the Hamilto- 
nian, is periodic in $ with the period of 2tt since the 
absolute phases of the basis wave functions used in the 
calculations have no effect on the diagonalization. How- 
ever, if we were to change the rotation phase adiabatically 
so that $ changes by 2tt, the wave function of an energy 
cigenstate could change, and along with it measurable 
properties of the system. The situation is analogous to 
what happens to a function with a branch cut starting 
from the origin of the complex plane: moving the com- 
plex argument of the function a full circle about the origin 
could put the value of the function to another branch. 

For instance, according to (36) and (37) one could ef- 
fectively flip the sign of J in the Hamiltonian by flipping 
the sign of every other basis function without any phys- 
ical consequences, but only if the number of the lattice 



sites is even. Otherwise the phase adjustments in the 
Hamiltonian over the lattice would not add up to an in- 
teger multiple of 2-k. This tells us that the physics may 
depend on the parity of the number of lattice sites. On 
the other hand, one could effectively flip the sign of J 
and adjust the energy spectrum accordingly by spinning 
the lattice up and thereby by adding -k to the rotation 
phase, but it is not obvious without further investigation 
what happens to the observable properties of the lattice. 



III. ONE-ATOM STATES 

We begin our investigation of the physics due to the 
rotation with the case of one atom in the lattice. If it 
were possible to extinguish atom-atom interactions, for 
instance with the aid of a Fcshbach resonance, many 
bosons may be put in the same one-particle state without 
any side effects. Such a many-atom sample could aid in 
the observation of the phenomena we will discuss. 

In the case of one particle the MTRF Hamiltonian (34) 
is diagonalized trivially. Each lattice momentum eigen- 
state annihilated by the boson operator B q is also an 
cigenstatc of energy with the characteristic frequency 



Now, given the RF momentum operator as shown in the 
Hamiltonian (11), we have 



-Jcos(q 



(38) 



The corresponding energy eigenstates are given by | ^ q ) = 
B^ |0), or, in the first-quantized representation in terms 
of the Wannicr-like function W, 

%(x) = *(<?, 0; *) = -Lj2 e^-^ n W(x - x n ) . (39) 



Let us now consider measurement of lattice momentum 
in a ring lattice by a variation of the methods that were 
successfully applied to measure lattice momenta in usual 
linear lattices [6, 24]. Imagine first that the ring lattice 
were straightened out, somehow in such a way that the 
boundary conditions etc. still work out as in the closed 
ring. The wave function of the state ^ q {x) is given by 
Eq. (39) in the MTRF, so we undo the momentum trans- 
lation and in the RF find 

%{x) = ~^Y. e l<in [e l ^- x ^/ a W{x - x n )] . (40) 



The wave function inside the square brackets is localized 
to a fraction of the lattice spacing a, and as such should 
have little effect on lattice momenta over the scale of the 
first Brillouin zone. Specifically, assume that the lattice is 
turned off so slowly that the details of the wave function 
on the scale of lattice spacing a get smoothed out, but 
fast enough that the structure over the length scale of 
the entire lattice L survives [6, 24]. With n = x/a, the 
wave function then becomes 



hd_ 
i dx 



mv I tjjq 



hq 



mv I ip q . 



(42) 



ip q oc e 



iqx/a 



(41) 



which identifies the measured kinetic momentum, mass 
times velocity, in the moving frame. Finally, applying the 
Galilean transformation, we see that in the LF the mo- 
mentum measured in this way would be Hq/a, simply the 
lattice momentum q translated into units of momentum. 
Therefore, as in the case with an ordinary stationary lin- 
ear lattice, the lattice momentum q may be converted to 
linear momentum and measured. 

In practice an operation that just straightens out a ring 
lattice appears implausible. Instead, think of a section 
of the ring that looks approximately like a straight piece. 
The atom, when released from different sections would 
have different directions of momenta, but still the same 
speed and energy The latter could then be measured. 
However, as design of experiments is not our aim, we 
will not pursue a more quantitative analysis. 

We have considered three different frames in our 
derivation: stationary lab frame LF, rotating frame RF, 
and momentum translated rotated frame MTRF. It may 
come as a surprise that, after we undo the momentum 
translation and are still formally in the RF, the lattice 
momentum q is the lattice momentum that would be 
observed in the LF. Nevertheless, an inspection of the 
Hamiltonian (11) shows why this is consistent: To get to 
the RF, we subtract the momentum due to the motion of 
the frame, mv, from the quantum mechanical momentum 
operator p = 5 -§- , which means that p should represent 
the momentum in the stationary LF. All of this is consis- 
tent with the observations surrounding Eqs. (15) -(17). 

On the other hand, when we quote values of energies, 
they are always in a moving frame, RF or MTRF. Mo- 
mentum translation has no effect on energy, as it comes 
with canceling transformations of both the state and the 
energy operator (Hamiltonian). 

We plot in Fig. 1 the characteristic frequencies of all 
energy eigenstates as a function of the rotation phase <f> 
for a lattice with N = 3 and 4 sites, upper and lower 
panels. There is a qualitative difference in the spectra 
in accordance with the observation that the parity of the 
number of states may make a difference. Since the plot 
presents ui q /J, for J > the energy increases from bot- 
tom to top on the vertical axes. 

Suppose first that, in fact, J < holds true, so that 
lowest energies arc at the top. Without rotation the 
ground state of a lattice with an odd number of sites 
is doubly degenerate. Moreover, the degenerate states 
have q ^ 0, i.e., present a flow along the lattice. If we 
prepare the ground state and measure the lattice momen- 
tum, we find one of two nonzero values at random. If the 
lattice has an even number of sites, the ground state is 
nondegenerate. It has the lattice momentum tt, which is 
physically equivalent to — w. Depending on the measure- 
ment scheme the result may be n or — n, but they count 
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FIG. 1. Energies of the lattice eigenstates for TV = 3 (upper 
panel) and N = 4 (lower panel) sites as a function of the site- 
to-site rotation phase cj>. The graphs are for lattice momenta 
q = (solid black line), q — 2-k/N (red dotted line), q — A-k/N 
(green dashed line) and q — 6n/N (blue dash-dotted line). 



as the same. 

On the other hand, for J > the ground state is non- 
degenerate and has zero lattice momentum, no matter if 
the number of lattice sites is even or odd. The sign of the 
hopping matrix elements may have dramatic qualitative 
consequences if the number of lattice sites is odd. 

Incidentally, in the examples we have studied analyt- 
ically and numerically, we always found J > 0. There 
may be a general principle dictating that for the kind of 
passive tunneling we have considered J > should hold 
true. But if so, active schemes that control the tunneling 
externally, say, using light-induced Raman transitions, 
are clearly not restricted to J > 0. 

As we have already noted, the question of the ground 
state is precarious once the lattice rotates. Let us posit 
that preparation to a very low temperature produces the 
lowest-energy state in the rotating frame. We also as- 
sume J > for the sake of the argument, and take 
the example with N = 4. Consultation of the lower 
panel in Fig. 1 then shows that, by varying the phase 
<j) over a 27r interval, for each q a region of <f> emerges 
such that the corresponding state q is the ground state 
and gets prepared. For instance, q = is a ground state 
in the interval 4>/2ir € [—1/8, 1/8], q = tt/2 in the inter- 
val 4>/2-k £ [1/8,3/8], and so on. Such switching of the 
nature of the ground state in a rotating ring lattice as 
a function of the rotation speed has been noted before, 
and in the many-atom case with weak atom-atom inter- 



actions it may lead to a Schrodinger cat superposition of 
macroscopic flow states [25, 26]. 

Thermal preparation presents an example about the 
roles of the rotation phases and system observables. Con- 
sider the ground state first around </> = 0, which in the 
MTRF reads 



1 



* (x) = -=52e- i + n W(x-x n ) 



(43) 



On the other hand, the ground state around </> = ir/2 has 
q = 7r/2, and the wave function can be written, among 
others, in the form 



*„ /2 (x) = -!= Y, e-«*-*M"W[x - x n ) 



(44) 



This is the same function of x as (43), albeit with the 
replacement cf> — > cf> — tt/2. As a function of <f> and x, the 
ground state for all <fi is the same as the ground state in 
the interval (t>/2n g [—1/8, 1/8) repeated in <j> with the 
period 2ir/N = tt/2. Nevertheless, even though in the 
MTRF the wave function repeats with the period A<j> = 
2tt/N, this is not necessarily the case with the observed 
quantities. In fact, lattice momentum as measured in the 
LF steps by 2tt/N every time the q label of the ground 
state switches. This would probably be the educated 
guess of most colleagues familiar with vortices in trapped 
gases. 

Things are different in an interesting way if one con- 
siders "slow" time dependent variation of the rotation 
velocity. We call such variation adiabatic. At this stage 
there are no interactions between the energy eigenstates 
and they (Bloch states) belong to different discrete sym- 
metries, namely, eigenvalues of the lattice translation op- 
erator U. We therefore assume that they are extremely 
robust even if a degeneracy is crossed when the rotation 
speed is varied in time. Put differently, when two states 
cross, they have a difference of an integer multiple of 2ir 
in their end-to-end phase winding. A continuous dynam- 
ics cannot discontinuously change the phase winding, so 
that transitions between the states are not possible. In 
Fig. 1 one would follow a curve of a given color with vary- 
ing <j>, and the lattice momentum observed in the LF is 
independent of the rotation speed of the lattice. 

It might seem peculiar that in a lattice an energy eigen- 
state such as the ground state cannot be wound up adi- 
abatically, but this appears to be a manifestation of the 
same phase rigidity that sustains a persistent current in 
a ring. We argued a while ago [27] that one has to cut 
the supcrfluid inside an essentially one-dimensional ring 
and thereby severe the continuity of the phase if one is 
to alter the state of circulation. The lattice potential 
cuts the ring, but not all the way. If the rate of change 
of the rotation phase $ is sufficiently small compared to 
the tunneling matrix element J, the "fluid" of the single 
atom can still adjust and respond as if it were continuous. 

Continuing from the preceding ground-state example, 
one could first prepare the ground state ^ v /2(x) in a 



rotating lattice corresponding to the rotation phase <f> = 
7r/4 and then wind down the rotation adiabatically while 
the state \E f 7r /2(^) stays this way. This is a way to prepare 
in principle any eigenstate of the stationary lattice. 

As one more item of the phenomenology of the rotat- 
ing lattice, let us discuss the fate of a wave packet. We 
use group velocity as the tool. The concept of group 
velocity is increasingly useful, the larger the number of 
lattice sites and the broader the wave packet (so that the 
role of dispersion is diminished) . We will not carry out 
the associated quantitative analysis for finite-size lattices 
or wave packets, but simply assume that group velocity 
yields useful qualitative predictions. It is 



dk 



a J sm(q — 



(45) 



where, in order to permit a direct comparison with the 
rotation velocity, we have included the lattice constant 
a explicitly, as in k = q/a. Group velocity governs the 
evolution of the spatial envelope of the wave packet, and 
is given by the expression (45) in both the RF and the 
MTRF. 

Take a localized stationary wave packet with q ~ 
prepared when the lattice does not rotate, and assume 
that the lattice is then set in motion adiabatically in such 
a way that the decomposition of the wave packet into its 
component q states is not perturbed by the process. In 
the rotating frames the wave packet therefore picks up 
the group velocity 



v g = —aJ sine 



1 ma 2 J 



■v, 



(46) 



where the latter form applies for a small rotation velocity. 
The group velocity is small compared to the rotation ve- 
locity, both because of the factor I/TV and because of the 
second factor that is essentially the ratio of the tunnel- 
ing frequency to the photon recoil frequency associated 
with an atom in the optical lattice. This means that the 
wave packet starts moving approximately, but not ex- 
actly, with the rotating lattice. Depending on the sign of 
the transition matrix element J, the wave packet either 
slightly lags or may even lead the lattice. We empha- 
size the curious contrast: lattice momenta are unchanged 
when the lattice is spun up, but the wave packet basically 
tracks the moving lattice. 

We have described several peculiar phenomena that 
should occur with a single particle or noninteracting par- 
ticles in a rotating ring lattice. Some of them, such as 
the scheme to prepare an arbitrary eigenstate of energy, 
were based on nontrivial assumptions. But then, we can 
turn the tables and say that an experiment would test 
the validity of these assumptions. 



IV. LATTICE DIMER 



in detail before [14], but some reorientation is in order 
here. While our emphasis was on the limit of an infinitely 
long lattice and the periodic boundary conditions were a 
matter of convenience, in a laboratory ring lattice the 
topology of the ring is, and a small number of sites could 
be, a physical reality, and may necessitate a numerical 
solution of the system. For the most part, however, our 
main emphasis is on the rotation phases, and when pos- 
sible we piggyback on our earlier analysis [14] of lattice 
dimers. 

Thus, we write the most general state of two bosons in 
the lattice in the lattice momentum representation as 



= E^)4 *L J°> 



P+q ~P-q 



(47) 



where |0) is the vacuum with no atoms present. Here P 
is the total lattice momentum of the dimer of sorts, the 
value of the conserved quantity 



P = T,qBlB q 



(48) 



We let the value P of the operator P range from — 2-zr 
to 2tt, so that ^P runs over the usual interval [— tt,tt) of 
lattice momenta. Also, \P need not be a valid lattice 
momentum, but ^P ± q have to be. This means that 
in the sum the lattice-momentum like quantity q may 
either run over legal "integer" lattice momenta, or it may 
be displaced from legal lattice momenta by a half-step 
ir/N. This is what we earlier termed "half-integer" lattice 
momenta. Either way, the sum over q in Eq. (47) runs 
over a set of N values so that \P + q and |P — q both 
run once over all permissible lattice momenta, with no 
two values separated by 2n or more. 

The coefficients A(q) govern the internal structure of 
the dimer. They remain to be determined. By the 
boson symmetry they can be, and are, chosen so that 
A(q) = A(—q). For convenience we regard the coeffi- 
cients A(q) as periodic over 2n. The states associated 

with the coefficients A(q) and A(-q), B\ b\ |0) 

2 p M 2 PTq 
are the same. As a result, the usual inner product of the 
states of the lattice system is expressed in terms of the 
expansion coefficients as 



ty A ,rl> B )=2Y d A*{q)B{q). 



(49) 



An explicit solution to the problem of a few-site lattice 
has to be tailored to the even/odd number of lattice sites 
and integer/half-integer lattice momenta, with different 
sets of possible values of the lattice momenta q [14], but 
we will not embark on an enumeration of the various 



We next proceed to the case of two interacting atoms 
in the lattice. We have discussed an analogous situation 



We now turn to the specifics of a rotating ring lattice. 
In the MTRF the time independent Schrddingcr equation 



gives an equation for the coefficients A{q) 

-2Jcos^P-^co S qA(q) + jJ2 A ^ = J A (iy 

(50) 
Defining the overall frequency scale Q(P, (j>) as 

Q(P,4>) = 2Jcos(±P-0), (51) 

and dimensionless variables representing energy and the 
strength of atom-atom interactions 



- E r- U 
the energy eigenvalue problem may be written 



1 E- 



1 



cos q 



1 
K 



(52) 



(53) 



The corresponding unit-normalized energy eigenstatcs 
are defined by the coefficients 



A(ui,q) = 



i + cos 5 



,C(w) 



E 

. 9 



(w + cosg) 2 



-1/2 



(54) 

We [14-16] and others [6-8, 12] have discussed the na- 
ture of the solutions before. There is a band of continuum 
states restricted to the interval of characteristic frequen- 
cies (— [Oj , |0|), and one bound state of the lattice dimcr 
that peels off from the continuum as the strength of atom- 
atom interactions is increased. The designations such as 
"continuum" are obviously only qualitative in the case of 
a finite and possibly even a small number of lattice sites. 
In the limit of a large number of sites, the the unsealed 
characteristic frequency of the bound state is 






t/ 2 



(55) 



As has also been noted many times before, the center- 
of-mass motion encompassed into the lattice momentum 
P does not completely separate from the internal degree 
of freedom, the variable q. The new feature here is that 
the rotation phase similarly, indirectly, influences the in- 
ternal structure of the lattice dimer. In view of Eq. (51), 
for two identical bosons the roles of center-of-mass lattice 
momentum and rotation phase may be interchanged as 
it comes to energetics. Any variation with respect to one 
may just as well be realized by varying the other. 

For attractive atom-atom interactions the bound 
state is the lowest-energy state, and by our statistical- 
mechanics assumption it gets prepared at zero temper- 
ature. Given the rotation phase </>, the lowest energy 
occurs for P ~ 2<j>, so that thermal equilibration at low 
temperature and with fixed rotation phase <f> will in gen- 
eral prepare a finite flow velocity for the molecule. 

Also, in the presence of atom-atom interactions the 
bound state is always separated from the continuum by 



a nonzero amount, so that adiabatic manipulations of the 
bound state are possible. One could prepare thermally a 
stationary lattice dimcr (P = 0) and then add adiabati- 
cally an arbitrary phase <f>. As far as the energetics of the 
dimer is concerned, this is physically equivalent to gen- 
crating the center-of-mass lattice momentum P = —2<p. 
This applies to the spatial structure of the dimcr. The 
rms size of the bound dimcr is An = ^/2\fl/U\ 2 [14], so 
that the dimer could be shrunk to a single site by choos- 
ing a rotation speed such that O = 0. Similarly, the 
spectroscopy of the dimer [14] may be controlled. 

The final item to understand is measurements of lattice 
momentum. As we have noted already, the label q even 
in the MTRF directly corresponds to a lattice momen- 
tum measured in the LF. Completely analogously to the 
earlier analysis of the stationary lattice [14], the prob- 
ability to find a lattice momentum q is proportional to 
\A(^P — q)\ 2 . The rotation phase does not directly enter 
this expression. In the limit N 3> 1 the unit-normalized 
(in the sense of the integral over q) probability density 
for the lattice momentum q is found to be 



fid) 



|K(P)P 



2lly/l + !C(Py 



{cos(q - \P) + sgn[/C(P)] y/1 +K(PfY 



with 



fC(P) = 



u 



2Jcos(iP 



(56) 



(57) 



Consider past experiments in which repulsively bound 
lattice dimers were produced with P ~ and the lattice 
momenta of the atoms were subsequently measured. The 
lattice momenta were predominantly found at the edges 
of the first Brillouin zone, q ~ 7r [6]. Suppose now that, 
after the repulsively bound pairs have been created, it 
would be possible to spin up the lattice adiabatically to 
the rotation phase (f> — n. By Eqs. (56) and (57), the 
effect is the same as reversing the sign of the atom-atom 
interaction. Correspondingly, the lattice momenta would 
then be found predominantly at the center of the Bril- 
louin zone. This should be contrasted to the observation 
that adiabatic variation of the rotation speed has no ef- 
fect on the lattice momentum of a single atom. 

A repulsively bound state (U > 0) lies above the con- 
tinuum band and an attractively bound state (U < 0) 
below. We doubt if any adiabatic method exists that 
converts one to the other, as in the process the state 
should move intact across the dissociation continuum. In 
our thought experiment, according to Eq. (55), the repul- 
sively bound state remains the highest-energy state the 
entire time while the lattice is set in motion. However, 
its signature changes to that of an attractively bound 
state. This is an indirect effect of the energetics. It is as 
if the rotation effectively changed the sign of the tunnel- 
ing matrix elements, and here it is the relative sign of the 
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tunneling matrix clement and the atom-atom interaction 
strength that counts. 

Still further opportunities for control open up if the 
atoms are not identical [9, 16], and in particular, if their 
masses are different. This implies that for the same ro- 
tation speed the rotation phases of the two atoms are 
different. We have discussed the Hamiltonian and the 
results that emerge for dissimilar species in detail before 
for a nonrotating lattice [16]. Here we briefly comment 
on the effects of the rotation phases, call them fa and 
fa for the two species 1 and 2. These could be bosons, 
fermions, or one of each. 

This time around the expansion coefficient A(q) de- 
scribing atom 1 with lattice momentum ^P + q and atom 
2 with lattice momentum ^P — q is no longer constrained 
to be an even function of q, and the correct form of the 
inner product now is as one might expect in the first 
place. 



(1&A,Vb)=X>*(<z)B(«). 



(58) 



With the definitions 



n = yjji + 

(3 = arctan 
the Schrodir 


Jl + 2Ji J 2 cos(P -fa-fa) 
\Ji-J2,_(P-fa-fa\ 


- \{fa - 


-fa), 

(59) 
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E 

u = — 

nn 

lger equation becomes 
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, + 0)4(g) + §y>(g')=u 


A(q). 


(60) 



matter only modulo the spacing between the states q, 
or modulo 2tt/N. On the other hand, a bound state 
appears in this system just as for identical bosons and 
likewise can be modified adiabatically. In such a case an 
arbitrary value of /3 may have physical relevance. 

As an example, let us consider lattice momenta. Just 
like in Rcf. [16], we may find in the limit JV»1 the prob- 
ability distribution, normalized to unity, that a measure- 
ment would find either atom 1 or atom 2 with the lattice 
momentum q, 



/1.2(g) 



|/C(P)p 



2^^i + jc{py 



K{P) = 



{sgn[JC(P)] v/l + /C(P) 2 + cos[ g - ±P ± 0\}* ' 
Via 



2J cos[UP - fa - fa)} 



(63) 



To limit the scope of the present exercise, we assume that 
the tunneling matrix elements are the same whereupon 
we have (3 = \{fa — fa), that it is possible to hold P + 
fa + fa = while varying fa and fa, and that K, is 
negative. Practicalities aside, this is not as far-fetched as 
it might sound, since the variable fa^ are only defined 
modulo 2-7T and their sum and difference may in general 
be varied independently even if they are proportional to 
one another. The species 1 and 2 then predominantly 
come out with the lattice momenta —fa and —fa. The 
novelty here is that until now we have not produced an 
example in which the position of the maximum of the 
lattice momentum distribution is somewhere else than at 
or ±7r. 



The branch of the explicit arctan function must be cho- 
sen judiciously so that the expression of /? is a continu- 
ous function of its variables. U12 is the strength of the 
interspecies interaction. The eigenvalue equation and its 
solutions may eventually be written 



N ^ 
1 

A(ui,q) ex 



1 



1 
\T ^ u + cos{q + (3) ~ K 



1 



us + cos(q + f3) 



(61) 
(62) 



A comparison with the case of two identical bosons 
first shows that, completely analogously, the sum of the 
rotation phases of the two atoms always gets subtracted 
from the center-of-mass lattice momentum. In the con- 
trol of the system, fa + fa and P are for the most part 
equally good knobs to turn. 

Next consider the eigenvalue problem (61), (62), for 
the sake of simplicity ignoring the dependence of K, on 
the rotation phases. The values of the parameter j3 then 



V. CONCLUDING REMARKS 

We have developed a Hubbard model for atoms in a 
rotating ring lattice carefully and in great detail in order 
to correctly identify the physical obscrvables. It turns out 
that thermal preparation and adiabatic variation of the 
rotation speed can be used to control the states of both a 
single atom and a dimcr of two atoms, occasionally with 
unexpected results. 

Our Hubbard model could be a rather cursory approx- 
imation under many experimental circumstances. For 
instance, it might happen that the motion of the atoms 
is not strictly confined to one spatial dimension, or a 
high rotation speed or atom-atom interactions render the 
tight-binding approximation questionable. However, the 
phenomena we have described arise from general princi- 
ples such as the difficulty of abruptly altering the phase 
winding of a quantum state around a ring, and should 
survive modest experimental imperfections. 
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